Photoproduction of Neutral Pions off Protons 
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Photoproduction of neutral pions has been studied with the CBELSA/TAPS detector in the 
reaction 7p —J- pTr" for photon energies between 0.85 and 2.50 GeV. The vr" mesons are observed in 
their dominant neutral decay mode: 7r° -^ 77. For the first time, the differential cross sections cover 
the very forward region, ^c.m. < 60°. A partial- wave analysis of these data within the Bonn-Gatchina 
framework observes the high-mass resonances Gi7(2190), Di3(2080), and Di5(2070). 

PACS numbers: 13.60.Le Meson production, 13.60.-r Photon and charged-lcpton interactions with hadrons, 
13.75.Gx Pion-baryon interactions, 14.40.Aq tt, K, and rj mesons, 25.20.Lj Photoproduction reactions 



I. INTRODUCTION 



Spectroscopy has long provided essential information 
for scientists trying to understand the nature of parti- 
cles. The spectroscopy of the hydrogen atom revealed 
the quantum nature of the interactions between particles 
at small distance scales and led to the development of 
quantum electrodynamics (QED). Exploration of the ex- 
cited states of hadrons opens a door into the nature of 
nonperturbative quantum chromodynamics (QCD). By 
mapping out the excitation spectra of hadrons, informa- 
tion is gained about their effective degrees of freedom and 
the effective forces that lead to confinement. A both par- 
ticularly challenging and interesting task is to understand 
the spectra of baryons. Although symmetric quark mod- 
els, describing the baryon with three constituent quark 
degrees of freedom, are good at predicting many prop- 
erties for many of the low-lying baryons (especially the 
ground states), there is a large disagreement between the 
number of excited states predicted by these models and 
the number that has been observed at and, in particular, 
above masses of about 1.8 GeV/c^. Attempts at finding 
a physical mechanism that can account for this discrep- 
ancy have not thus far yielded a satisfactory answer, and 
the higher mass regions of the excitation spectra remain 
largely unexplored: Most known baryon resonances have 
masses smaller than 2 GeV/c^ and were discovered in 
elastic ttA'^ scattering experiments. Many of the missing 



or hitherto unobserved baryon resonances are predicted 
to couple only weakly to the ttA^ channel but exhibit large 
couplings to ^N . In recent years, many different labo- 
ratories worldwide (Jefferson Laboratory, ELSA, MAMI, 
GRAAL, SPring-8, etc.) have measured differential cross 
sections and polarization observables for a large variety 
of final states in electromagnetically induced reactions. 
A good review is given in [l|. 

Although many of the missing baryon resonances may 
have a small coupling to ttTV, it is still important to study 
pion photoproduction. Doing so may confirm or reject 
resonances seen in elastic pion-nucleon scattering, and 
it could still provide evidence for unexplored resonances. 
New resonances found in reactions such as 7iV -^ ttN are 
expected to have masses larger than about 1.8 GeV/c^, 
although the higher-mass resonance contributions arc ex- 
pected to be more important in double-meson photopro- 
duction, by sequentially decaying via emission of a tt or 
77 meson and populating intermediate states. 

Excited baryons are broad and overlapping, making it 
difficult to extract their properties by studying isolated 
reaction channels and cross sections alone. Polarization 
observables can provide access to small resonance contri- 
butions to a final state, so that missing resonances may 
still be found by studying these in pion photoproduction. 
However, both unpolarized cross sections and polariza- 
tion observables are needed to maximally constrain and 
uniquely determine the scattering amplitude for a given 
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FIG. 1: (Color online) Cross section of the CBELSA/TAPS setup (not to scale). The electron beam delivered by the accelerator 
ELS A enters from the left side and hits the copper radiator (for the data presented here) in front of the tagger magnet. 



process. Additional physical constraints arc provided by 
performing an analysis of these observables for multiple 
reaction channels simultaneously, and to succeed in its 
description of the contributing resonances, the combined 
analysis of these data must accurately describe the non- 
resonant behavior of the observables. Previous experi- 
ments have shown that t-channel contributions are im- 
portant in pion photoproduction at high energies, where 
the missing resonances are expected. Thus further study 
of the pion photoproduction cross section, particularly at 
forward angles, is very interesting. 

Recent measurements of differential cross sections in 
pion photoproduction have made significant contribu- 
tions to the world database, but these measurements have 
limited coverage at very forward center-of-mass polar an- 
gles of the pion, 6lc.m.- The recent CB-ELSA g, Q and 
CLAS data Q on 7p — )► pir^ cover angles down to approx- 
imately 30° and 40° in 0c. m., respectively. They show dis- 
crepancies, where the respective differential cross sections 
are observed to increase at different rates as one moves 
to more forward angles for the incoming photon energy 
range 1300 to 2000 MeV. Our new CBELSA/TAPS data, 
presented here, extend the coverage of the forward region 
and will help understand the behavior of the tt" photo- 
production cross section in this kinematic region. 

In this paper, we present differential cross sections for 
the reaction: 



7P — >■ pTT 



where 7r° — > 27. 



(1) 



We studied, in particular, the forward region; a similar 
study was published for the beam asymmetries S in the 
same reaction Q. 

The paper has the following structure. Section |ll] sum- 
marizes recent results on tt*^ photoproduction. A brief 
introduction to the CBELSA/TAPS experimental setup 
is given in Section lllll The data reconstruction and event 
selection are discussed in Section ITVl and the extraction 



of cross sections is described in Section [V] Experimental 
results are finally presented in Section I VII and results of 
a partial- wave analysis are discussed in Section FVlII 



II. PREVIOUS RESULTS 

Cross section data for tt" photoproduction were ob- 
tained and studied at many different laboratories over a 
wide kinematic range [2|-lj, l6l-[l8l| . A review of the main 
data sets pubhshed before 2005 and a corresponding com- 
parison of their coverage in energy and solid angle can be 
found in [3|. The new data presented in this paper ex- 
tend previous data by covering the very forward polar 
angles with an electromagnetic calorimeter, TAPS. 

In 2005, the CB-ELSA collaboration at Bonn, Ger- 
many, presented high-statistics results on the photopro- 
duction of tt" mesons using the Crystal Barrel detector 
for incident photon energies from 300 to 3000 MeV 01. 
Above energies of 1300 MeV, the data cover an angu- 
lar range from about 30° to 140°, which corresponds to 
about -0.75 < cos 6*0. m. < 0.85. 

The GRenoble Anneau Accelerateur Laser (GRAAL) 
collaboration at ESRF in Grenoble, France, measured 
differential cross sections over a wide angular range 
from 550 to 1500 MeV [H, though still limited to 
cos 0c. m. < 0.7 or approximately 0c. m. > 45°. At 
the GRAAL facility, Compton backscattering of low- 
energy photons off ultrarelativistic electrons reached al- 
most 100 % beam polarization at the Compton edge. 

More recently, the LEPS collaboration at SPring-8 in 
Hyogo, Japan, published results for higher photon ener- 
gies between E^ = 1500 and 2400 MeV 0. For the first 
time, the data cover the tt" backward angles between 
-1 < cos^c.m. < -0-6 or about 125° < 0c.m. < 180°. 
Backward-Compton scattering was applied using Ar-ion 
laser photons with a 351-nm wavelength. 
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FIG. 2: (Top) Schematic drawing of the liquid-hydrogen 
target, scintillating-fiber detector, Crystal Barrel and TAPS 
calorimeters. (Bottom) Front view of TAPS: The left side 
shows the logical segmentation for the LED-low trigger; the 
right side shows the logical segmentation for the LED-high 
trigger (see text for more details). 



At Jefferson Lab, Newport News, VA, differential cross 
sections for the reaction were measured using the CEBAF 
Large Acceptance Spectrometer (CLAS) and a tagged 
photon beam with energies from 0.675 to 2.875 GeV [J]. 
The cross-section data cover an angular range from 
^c.m. ~ 40° to 150°, corresponding to approximately 
-0.85 < cos6ic.m. <0.75. 



III. EXPERIMENTAL SETUP 

The experiment was conducted using the Crystal Bar- 
rel [13 and TAPS [20, l2l[ electromagnetic calorimeters 
at the electron stretcher accelerator ELSA [22|, which is 
located at the University of Bonn in Germany. Figure [T] 
contains a schematic of the experimental setup. 

A 3.175-GcV electron beam was extracted (by slow res- 
onant extraction) from ELSA. The electrons were passed 
through a thin copper radiator, which had a thickness 
of (3/1000) Xji (radiation length), where a fraction of 
them underwent bremsstrahlung. These electrons were 
deflected by a dipole magnet into the tagger detector sys- 
tem (tagger), which consisted of 480 scintillating fibers 
above 14 scintillation counters in a partly overlapping 
configuration. Good knowledge of the magnetic field and 



the position of the hit in the tagger allowed the mo- 
mentum and energy of each electron to be determined. 
The tagged-photon energy is given by the difference be- 
tween the incoming electron energy and the electron en- 
ergy after bremsstrahlung. Electrons that did not un- 
dergo bremsstrahlung were deflected at small angles into 
a beam dump located behind the tagger. Tagged photons 
had energies ranging from 0.5 to 2.9 GeV; the higher en- 
ergy part of this bremsstrahlung spectrum had an energy 
resolution of about 2 MeV, while the lower part had an 
energy resolution of approximately 25 MeV. 

The energy calibration of the photon tagging system 
was determined in simulations using the known positions 
of the scintillating fibers and the measured magnetic field 
map, and approximated by a polynomial. Measurements 
with a known (lower- energy) ELSA beam allowed for a 
cross-check of the calibration; this involved guiding a low- 
current beam directly into the tagger, while holding the 
magnetic field constant. Corrections to the original poly- 
nomial were determined from this procedure. 

Tagged photons continued along the beam line toward 
a liquid-hydrogen target, which was 5 cm in length and 
3 cm in diameter, located at the center of the Crystal 
Barrel (CB) calorimeter (see Fig. [T]). Some of these pho- 
tons interacted with protons within the target to pro- 
duce a proton and multiple photons in the final state, 
the photons resulting from the decay of neutral mesons. 
Final-state photons were detected by the two electro- 
magnetic calorimeters (CB and TAPS), which together 
almost covered the full 47r solid angle (> 98% coverage). 
The CB detector in its CBELSA/TAPS configuration of 
2002/2003 consisted of 1290 CsI(Tl) crystals, which had 
a trapezoidal shape and were oriented toward the target, 
with a thickness of 16 X^, and covered the polar angles 
from 30° to 168°, while TAPS consisted of 528 hexagonal 
BaF2 crystals with a length of approximately 12 Xji-, and 
covered the polar angles from 5° to 30° ; both covered the 
full azimuthal circle. TAPS was configured as a hexag- 
onal wall and served as the forward endcap of the CB. 
The CB originally covered polar angles down to 12° and 
was made of 26 rings [l^. Downstream rings 11 — 13 of 
the CB were removed to allow for a tight fit with TAPS 
(see top of Fig. [21). Forward-going protons were detected 
by plastic scintillators (5 mm thick) located in front of 
each TAPS module; the other protons were detected by a 
three-layer scintillating fiber detector, which closely sur- 
rounded the target [23] . 

Photons that did not interact with the target material 
continued further along the beam line through a hole in 
TAPS and finally arrived at the beam monitor, a total 
absorption Cerenkov counter, which consisted of an ar- 
ray of nine lead glass crystals. The number of photons 
counted here was vital to the determination of the photon 
flux. 

The first-level trigger was provided by the fast re- 
sponse of the TAPS modules. A cellular logic called FAst- 
Cluster Encoder (FACE), which counted the number of 
clusters in the barrel, formed the second-level trigger. 
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FIG. 3: (Color online) Invariant 77 mass spectra for the reaction "fp — >■ P77 for E.y £ [875, 900] MeV and for (a) 6c. m. G [0°, 10°], 
(b) 6c. m. £ [40°, 50°], and (c) Oc.m. £ [160°, 170°]. CL cuts were applied at 10~^. The red solid line denotes the signal Gaussian, 
the green dotted line the Landau distribution for the background description, and the blue dotted line shows the background 
subtracted distribution. The background is essentially negligible for 6c. m. > 50°. 



The trigger required at least one hit in TAPS above a 
high-energy threshold (LED-high) and at least one FACE 
cluster, or it required at least two hits in TAPS in differ- 
ent segments above a low-energy threshold (LED-low). 
The logical LED-low and LED-high segmentation of the 
trigger is shown in Fig. [2l Topologically, reaction ([1]) 
events satisfying the first condition have either a fast pro- 
ton in TAPS (and a backward-going pion in the CB) or 
one high-energy photon in TAPS and one photon in the 
CB (from a forward-going pion) , while those that satisfy 
the second condition have a pion in TAPS, together leav- 
ing a hole in the pion acceptance at intermediate angles. 



IV. PREPARATION OF FINAL STATE 

The data presented here were recorded from October 
2002 until November 2002 in two run periods at ELSA 
beam energies of 3.175 GeV. For this analysis, only the 
photon energy range up to 2.55 GeV was used since 
the tagger consisted of scintillating fibers covering en- 
ergies up to this energy and providing timing informa- 
tion, while simple hit/miss wires covered photon ener- 
gies above 2.55 GeV. These data were also used for pre- 
viously published analyses on a variety of different final 
states 0411]. 

This section discusses the event reconstruction and se- 
lection of the TT^ photoproduction channel: reaction ([IJ. 
A total number of 2.21 x 10^ tt" events was observed, 
covering invariant masses from 1570 to 2360 MeV/c^. 



A. Event Reconstruction 

The event selection was carried out in close analogy to 
the study of the 77 photoproduction channel. Further de- 
tails on the event reconstruction and the use of kinematic 
fitting in this analysis can be found in [26[. 

Events were analyzed for which only the two final-state 
photons (from the decay of the pion, 7r° — ;■ 77) or all 
final-state particles (77 -|~ p) were detected. The protons 
were identified by using the photon- veto BaF2 detectors 
of TAPS or the inner scintillating fiber detector inside the 



CB. A CB cluster is assigned to a charged particle if the 
trajectory from the target center to the barrel hit forms 
an angle of less than 20° with a trajectory from the target 
center to a hit in the inner detector. The resolution for 
proton clusters in the calorimeters is worse than that for 
photon clusters because, on average, proton clusters are 
much smaller than photon clusters. Proton identification 
was only used to separate protons from photons, but the 
proton four-vector was determined from the event kine- 
matics through missing-proton kinematic fitting. 

Time accidental background was reduced by demand- 
ing a prompt coincidence between an electron in the tag- 
ger and a particle in TAPS; if this particle was a photon, 
it required a time difference within ±3 ns, while for a pro- 
ton a time difference within -5 to -1-15 ns was required. 
Two to three tagger photons, at most, survived the initial 
timing cut. A missing-proton kinematic fit was used to 
separate the best initial-state photon from the remain- 
ing candidates in the event. This IC kinematic fit to 
the hypothesis 7p — >■ p miss 77, which only required four- 
momentum conservation, was performed for each initial 
photon candidate of an event passing the initial timing 
cut, and the candidate with the largest confidence level 
(CL) or x^ probability was selected. Finally, only events 
with a CL > 1 % were accepted for further analysis. No 
mass constraints were used in the kinematic fitting. 

In the final step of the event reconstruction, we re- 
moved events that are not p(7r° -> 77) final states (back- 
ground) from the event sample. This was done by inter- 
polating the shape of the background in the signal region 
based on the distribution of events outside the signal re- 
gion. The data were binned, and maximum likelihood 
fits of a signal plus background hypothesis to the data 
were performed. A Gaussian line shape was used for 
the signal, and a Landau distribution was used for the 
background. Typical 77 mass distributions are shown 
in Fig. El for E^ £ [875,900] MeV. The size of the back- 
ground at very forward angles is substantial [Fig.[3](a)]; it 
is much smaller, on average, in the backward bins at less 
than 4 % [Fig. [3](c)] . Mass resolutions (standard devia- 
tions of Gaussian signals) extracted from the data agree 
within 1 MeV/c^ with those from Monte Carlo simula- 
tions. 



B. Monte Carlo Simulations 

The performance of the detector was simulated in 
GEANT3-based Monte Carlo studies. We used a pro- 
gram package that has been built upon a program de- 
veloped for the CB-ELSA experiment. The Monte Carlo 
program accurately reproduces the response of the TAPS 
and CB crystals to photons. 

The acceptance for reaction ([T]) was determined by sim- 
ulating events that were evenly distributed over the avail- 
able phase space. The tagging and timing of initial-state 
photons were not simulated. The Monte Carlo events 
were analyzed by using the same reconstruction criteria 
as were also applied to the (real) measured data. The 
same IC-hypothesis was tested in the kinematic fits, and 
events were selected with the same CL cuts. The ac- 
ceptance is defined as the ratio of the number of recon- 
structed to generated Monte Carlo events, 



7P -^ p X 



N, 



■cc,MC 



N, 



gcn,MC 



{X^n"). 



(2) 



In the analysis presented here, we have applied an ac- 
ceptance cut of at least 5% on {E^, cos 0c. m.) bins, as 
well as an acceptance cut of at least 1 % on {E^, dan. ) 
bins, and removed the other data points from the analy- 
sis. Different values of the acceptance cut were used for 
the two binnings to provide comparable statistics over a 
large angular range. 



6.022-10^^ mol~i is the Avogadro number, and L = 5.275 
cm is the length of the target cell. The factor of two ac- 
counts for the molecular composition of hydrogen (H2). 

The number of events in an (£^^, cos^c.m.) or {Ej, 
^c.m.) bin comprises events with two or three final-state 
particles (at least 27's). The proton can be "miss- 
ing" , but events with and without a detected proton are 
treated in the same way in the event reconstruction. For 
cos 0c. m. < 0.0, the event kinematics requires that the 
proton is used in the (TAPS) trigger. Thus, the detec- 
tion efficiency of forward-going protons needs to be rea- 
sonably well understood. This has been shown earlier 
for CBELSA/TAPS results on 77 photoproduction; more 
details on the proton trigger can be found in [2g . 

The interval of the solid angle is given by A 51 = 
27rA cos(6'c.m.), with A 008(00. m.) describing the width of 
the angular bins. It was chosen to be 0.1 for the 7r° data 
presented here. For the representation of the differential 
cross sections in terms of the angle 0c. m. directly, a bin 
width of 10° was chosen. Energy bins were defined by 
considering statistics and ensuring a good comparability 
with other experiments. A total of 42 bins is presented 
in energy steps of 25 McV for E^ G [850, 1300] MeV 
and 50 McV for E^ G [1300, 2500] MeV. 

The number of observed tt'^ mesons needs to be cor- 
rected for unseen decay modes. The partial-decay 
branching fraction used to correct the measured cross 
sections is taken from [13 : BR(7r° -^ 77) = 0.98823 ± 
0.00034. 



V. DETERMINATION OF CROSS SECTIONS 



Normalization 



The differential cross sections for this analysis have 
been determined according to 



do- 
dH 



N^o_ 



1 



1 



1 



A^o^^^ N^pt An EjiI 



(3) 



r total 



where 



N, 



A, 



'->77 

AD. 



target area density 

number of reconstructed data events 



{E, 



or (S^,6lc.m.) bin 



bin 



number of beam photons in an (E-y) bin 
acceptance in an {E.y,cos6c.m.) or {Ej,9c.i 
solid-angle interval AQ, — 27rAcos {9 cm.) 
decay branching fraction. 



The target area density, that is, the number of atoms in 
the target material per cross-sectional area (orthogonal 
to the photon beam), is given by 



Pt =2 



p{R2)NaL 

Mmol(H2) 



= 2.231 • lO^Vb" 



(4) 



where p(H2) — 0.0708 g/cm'^ is the density, M^^i 
2.01588 g/mol is the molar mass of liquid H2, Na 



The CBELSA/TAPS tagging hodoscope consisted of 
480 scintillating fibers above 14 partially overlapping 
scintillation counters (tagger bars). The photon flux was 
measured directly in the experiment and determined ac- 
cording to 



N^ 



a-P-y-N 



fiber 
scaler ' 



(5) 



where the parameter a accounts for the (fiber)-cluster 
reconstruction in the tagger, which has to be performed 
in the same way as for real hadronic events. Pj, the pho- 
ton definition probability, denotes the probability that a 
real photon is emitted along the beam axis in the tag- 
ger and traverses the liquid hydrogen target; it is deter- 
mined from Tagger- Or- Runs, separate data runs utiliz- 
ing a minimum-bias trigger. A'^gc'aicr ^^ ^^^^ number of 
free hardware counts for the individual fibers corrected 
for the lifetime of the detector. Scalers are recorded in 
scaler events, which were accumulated with a minimum- 
bias trigger at a rate of 1 Hz during regular data taking. 
This trigger required only a hit in the tagger and was 
thus independent of hadronic cross sections. 

The Pj error is assumed to dominate the total error of 
the photon flux, which has a strong dependence on the 
efficiency of the beam monitor at the end of the beam 
line (Fig. [T]). By varying the background subtraction 



of noncoincidcnt taggcr-bcam monitor hits, Py was de- 



termined to be 0.639 ± 0.002. 



± 0.05 = 



This 



' stat. -^ yj.\jtj sys. 

value agrees with determinations from muhiplc Tagger- 
Or-Runs at different incoming photon rates. An overaU 
error of 10% has been assigned to the photon flux de- 
termination and has been included in Figs. |4] and O 



B. Systematic Uncertainties 



The statistical errors are determined from the number 
of events in each {E^, cos dan.) or (Ey,6c.m.) bin. Sta- 
tistical errors are shown for all data points; systematic 
uncertainties are given as error bands at the bottom of 
each plot. 

The systematic errors include uncertainties from the 
position of the liquid hydrogen target and a possible 
offset of the photon beam. The position of the target 
cell was determined from kinematic fitting by comparing 
the off-zero displacement of different pull distributions to 
Monte Carlo simulations. It was found to be shifted up- 
stream by 0.65 cm [31 . The corresponding systematic er- 
rors were determined by varying the target position in the 
Monte Carlo (±1.5 mm) and evaluating changes in the 
reextracted differential cross sections. The errors show 
an angular dependence, but are 3-4% on average and 
< 7% at most in the very forward region, 6'c.m. < 50°. 
The photon beam was assumed to be shifted by less than 
2 mm off axis at the target position. The errors of the 
decay branching fractions are negligible. 

The reconstruction of photons from the decay of neu- 
tral mesons and the identification of final states requires a 
sequence of cuts including the use of kinematic fitting. In 
[26l | , it was discussed in detail that the reconstruction of 
77 mesons from final states with two and six photons leads 
to compatible results for the 77 differential cross sections. 
The reconstruction of tt^ (— >■ 77) mesons is very similar 
to the reconstruction of ?/ mesons in the same two-photon 
decay mode and for the 7r° analysis discussed here, the 
same data sets, software packages, and calibration were 
used. This shows the good understanding of the detector 
response to multi-photon final states. An overall ± 5.7% 
error is assi gne d to the reconstruction efficiency as de- 
termined in |28| . An additional 3 % systematic error ac- 
counts for the slightly different effects of CL cuts on data 
and Monte Carlo events. The uncertainty of the proton 
trigger has been estimated from the small disagreement 
of the differential i] cross sections using the 7/ — >■ 27 and 
77 —5- 37r° — > 67 decay channels for E-y < 1 GcV and 
cos 0c. m. < 0.0 where the proton trigger dominates over 
triggers caused by photons. 

All the discussed contributions to the systematic error 
have been added quadratically to give the total system- 
atic uncertainty for this reaction. 



VI. EXPERIMENTAL RESULTS 

Figure |4] shows the jp — > p7r" differential cross sections 
plotted versus cos 6'c.m. of the 7r°, using cos 6'c.m. bins of 
width 0.1. To facilitate the comparison with our previous 
CB-ELSA data, we have chosen energy bins of 25 MeV for 
the energy range E^ G [850, 1300] MeV and energy bins 
of 50 MeV for the energy range E^ £ [1300, 2500] MeV. 
The data cover the very forward region. We have ex- 
cluded those data points in this analysis that show a 
Monte Carlo acceptance of < 5%. The error bars are 
statistical only; the total systematic error is given as 
error bands at the bottom of each plot. The angular 
distributions exhibit a region of zero acceptance around 
approximately —0.3 < cos 0c. m. < 0.3, depending on en- 
ergy, because of the trigger conditions during the data 
taking. We observed a small drop of the cross section 
for Ej € [1075, 1100] MeV, owing to an increased photon 
flux determined for one of the tagger fibers defining this 
energy bin; the reason for this increase is unknown, and 
thus these data are not shown. Resonance production is 
clearly observed up to high energies, as can be seen from 
variations in the differential cross sections over the full 
angular range up to high energies. Above E~^ > 1.5 GeV, 
the development of a forward peak indicates important 
contributions from i-channel exchange, which becomes 
dominant above incoming energies of 2.2 GeV. 

Figure [5] shows the same jp —J- p7r° differential cross 
sections plotted versus Oc.m. of the 7r° to illustrate the 
improvement of the data coverage in the forward direc- 
tion. The width of each 6'c.m. bin is 10°. This change 
of the data representation and the modified binning re- 
quired a new analysis and subtraction of the background 
under the tt^ peaks in the mass spectra as well as a new 
investigation of the systematic errors. Again, the total 
systematic error is given as error bands at the bottom of 
each plot. For this representation of the data, we have ex- 
cluded data points that show a Monte Carlo acceptance 
of < 1%. Figure [HI shows the differential cross sections at 
three different incoming photon energies to facilitate the 
comparison with a variety of other data sets. 



A. Comparison with CB-ELSA and CLAS Data 

A more detailed comparison of the results in the form 
of cross section ratios of this analysis to other data sets 
is shown in Fig. [T] The red open circles represent the 
forward direction, cos0^o > 0, and the solid blue squares 
the backward direction, cos6'^o < 0. In Fig. [7](a), we 
compare the new CBELSA/TAPS data of this analysis 
with the previous CLAS datajj] and in Fig. [7](b), with 
the previous CB-ELSA data [4]. In Fig. [7](c), a compar- 
ison between CB-ELSA and CLAS is shown. Since the 
earlier CB-ELSA data switch from a 50-MeV to a 100- 
MeV energy binning at high energies owing to statistical 
limitations, the distributions in panels (b) and (c) are cut 
off at an incoming photon energy of 2.3 GeV. 
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FIG. 4: (Color online) Differential cross sections for "fp — >■ pTr" from CBELSA/TAPS (•) plotted versus cosSc.m.- For compar- 
ison, CB-ELSA data Ql are represented by (■) and CLAS data [J| by (o). Note that the CLAS cross sections were extracted 
for 50-MeV-wide energy bins; for E~f < 1.3 GeV, the same CLAS results are thus shown in the corresponding 25-MeV-wide 
bins chosen for this analysis. The solid blue line shows the solution of the Bonn-Gatchina partial- wave analysis [23], and the 
dashed black line represents the current SAID solution (SP09) [SOJ]. The data points include statistical errors only; the total 
systematic error is given as error bands at the bottom of each plot. 
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FIG. 5: (Color online) Differential cross sections for 7p — ^ pn" from CBELSA/TAPS (•) plotted versus Sc.m.- For comparison, 
CB-ELSA data Q] are represented by (■) and CLAS data [j| by (o). Note that the CLAS cross sections were extracted for 
50-MeV-wide energy bins; for _E^ < 1.3 GeV, the same CLAS results are thus shown in the corresponding 25-MeV-wide bins 
chosen for this analysis. The solid blue line shows the solution of the Bonn-Gatchina partial- wave analysis [2^, and the dashed 
black line represents the current SAID solution (SP09) [30| . The data points include statistical errors only; the total systematic 
error is given as error bands at the bottom of each plot. 



* (a) 1350- 1400 MeV 

*4 



(b) 1750 -1800 MeV 




(c) 1900 -1950 MeV 



50 100 150 



50 100 



50 100 150 e 



cm. 



FIG. 6: (Color online) The diflterential cross sections for 7p -^ p-K° at E-, £ [1350, 1400] MeV (left), E^ G [1750, 1800] MeV 
(center), and E^ £ [1900, 1950] MeV (right). The angle shown is the pion center-of-mass scattering angle. The experimental 
data are from the CBELS A/TAPS measurements presented here (filled circles •), CB-ELSA 0] (open stars ii), GRAAL 
(left) [Tsj (open boxes D), SPring-8 (center and rig ht) [l6l ] (open boxes D), CLAS [3] (open circles o), Althoff et al. (left) 
[17| (open triangles /,), and Brefeld et al. (right) [Ig] (open triangles a). The dashed lines represent the current SAID solution 
(SP09) [3y], and the solid lines represent the latest solution of the Bonn-Gatchina PWA [23], which includes the new data 
presented here. The total systematic error is given again as an error band at the bottom of each figure. 



Figure [8] shows the excitation functions for the reac- 
tion 7p — >■ pn'^ at fixed values of the pion center-of- 
mass scattering angle using a logarithmic scale. We also 
added available data points in the backward direction 
from LEPS [l^. The agreement with LEPS is good and 
confirms the weak energy dependence of the cross sec- 
tion found above 2 GeV and at very backward angles. A 
very similar, mostly flat distribution is observed at very 
forward angles. 

In a direct comparison, our new data show an over- 
all good agreement with the older CB-ELSA data within 
systematic errors [Fig.[7](b)]. While our new data show a 
steeper slope of the forward (i-channel) rise than the ear- 
lier CB-ELSA results, the most forward CB-ELSA point 
at 9c.m. ~ 35° is almost perfectly matched for all en- 
ergies, indicating the very good agreement of the overall 
strength of the forward peak (Fig. [51). The few red points 
in Fig. [7|(b), which are systematically higher between 
1.5 and 2.3 GeV, indicate the different slope observed in 
this analysis. A less rapid increase of the differential cross 
section in this kinematic region, 0.3 < cos^jrO, which was 
implied by the recent CLAS data and suggested in [J], 
can thus not be confirmed in this analysis. 

Below E^ < 1.3 GeV, the overall discrepancies be- 
tween CBELSA/TAPS and CLAS are smaller than about 
20%, and thus, both data sets are consistent within sys- 
tematic errors [Fig. [7|(a)]. At higher energies, the for- 
ward direction agrees fairly well between the two data 
sets and appears to be consistent within systematic er- 
rors. Only a somewhat higher band of forward points 
indicates again the steeper slope of the forward peak ob- 
served in this analysis. However, agreement with CLAS 
for 0.6 < cos 6*^0 is fair to good (Fig. [8]). Note that addi- 
tional systematic normalization errors of the CLAS and 
CB-ELSA data (discussed in [3, |j] , but not included in 



the published data) are not considered in this compar- 
ison. Larger systematic differences up to 50 % at and 
above Ej « 1.5 GeV remain only in the backward direc- 
tion. This confirms the findings of the previous CB-ELSA 
analysis, which also observed a higher differential cross 
section in the backward region relative to the CLAS data 
[Fig. [71(c)] . Older data from Bonn [13] also support these 
higher cross sections [Fig. [n](a)]. 

This forward-backward disagreement suggests a 
non-trivial scale discrepancy between CLAS and 
CBELSA/TAPS, and we conclude that a normalization 
effect can be ruled out. Most of the recent data sets are 




FIG. 7: (Color online) Ratios of differential cross sections 
from this analysis to CLAS H (A) and to CB-ELSA [3 (B), as 
well as from CB-ELSA to CLAS (C). The red open circles rep- 
resent the forward direction (cos 6^o > 0) and the blue squares 
the backward direction (cosS.^o < 0). For E^ < 1.3 GeV, the 
same CLAS results have been used for the corresponding 25- 
MeV-wide bins chosen for this analysis. 
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FIG. 8: (Color online) Fixed-angle excitation functions for the reaction 7p — >■ ptt"; the angle shown is the pion center-of- 
mass scattering angle. The experimental data are from the new analysis presented here (filled circles •), CB-ELSA Ql (■), 
CLAS [J] (open circles o), and LEPS [lg| (filled -k). Note that the CLAS cross sections were extracted for 50- Me V- wide energy 
bins; for E-y < 1.3 GeV, the same CLAS results are shown in the corresponding 25-MeV-wide bins chosen for this analysis. 



not statistically limited, and a reason for the observed 
differences is not known. Since all the discussed discrep- 
ancies show the overall level of uncertainty for the reac- 
tion jp -^ pn'^ and this has a significant impact on the 
extraction of resonance couplings, this issue needs to be 
further addressed and resolved by future experiments. 



VII. 



PARTIAL- WAVE ANALYSIS 



A partial- wave analysis (PWA) has been performed 
within the Bonn-Gatchina PWA framework and includes 
the data presented here. For the discussion of contribut- 
ing resonances, we use the conventional naming scheme 
ofthePDG M: 



iV(mass of the resonance) L 21, 2.) 



(6) 



where / and J denote isospin and total spin of the res- 
onance, respectively, and L the orbital angular momen- 
tum in the decay into a nucleon and a pion. The scat- 
tering amplitudes in the PWA for the production and 
decay of baryon resonances are constructed in the frame- 
work of the spin-momentum operator expansion method 



suggested in |31l - l33| . The approach is relativistically in- 
variant and allows one to perform combined analyses 
of different reactions imposing analyticity and unitar- 
ity directly. The analysis method is further documented 
in [SJ, [33 . The experimental database used for the PWA 
is described in detail in [29| and includes the number of 
data points for each reaction as well as the weight for each 
data set used in the fits. It takes into account almost all 
important sets of photo- and pion-induced reactions, in- 
cluding three-body final states. The latter were included 
in event-based likelihood fits for the PWA solution dis- 
cussed in [2^. A full description of the experimental 
database goes beyond the scope of this paper, and thus 
we limit the discussion to tt'^ photoproduction data. 

For the jp —>■ pn^ differential cross sections, data for 
the invariant mass region above 1600 MeV/c^ are avail- 
able from GRAAL [l|, CB-ELSA 0,13, and CLAS 0. 
However, the GRAAL data are available only up to 
1950 MeV/c^ and the CLAS data cover only a rather 
restricted angular range at large masses. The new 
CBELSA/TAPS data cover the forward and backward 
angular region up to 2370 MeV/c^ and therefore pro- 
vide important information about states in the third and 
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fourth resonance regions. In particular, the new data 
are rather sensitive to contributions from high-spin reso- 
nances located in the mass region 2050 — 2150 MeV/c^. 
The behavior of the differential cross sections at very ex- 
treme angles should be a good indicator for the presence 
of such states. The sensitivity of the PWA solution to 
a particular resonance is studied by omitting or includ- 
ing it in the fits and evaluating the fit quality in terms 
of the resulting ^ values. The observation of a clear 
minimum in a mass scan gives further confidence about 
a contributing resonance. In such a scan, the mass of a 
resonance is chosen in several steps in the vicinity of the 
optimum value and fixed in the fit. The mass is then 
plotted versus the change of the x^ value from the best 
fit. Figure [9] shows a typical scan from this analysis. 

An overall scaling factor of 14% is needed in the 
PWA to achieve quantitative consistency of the data with 
the previous measurements of the CB-ELSA collabora- 
tion [3|. It is worth noting that a similar factor was 
needed to provide a compatibility between the CB-ELSA 
and CLAS data Q. In that case, however, the CLAS data 
needed to be scaled up by 5 %. Fitting only such a scale 
factor, we obtained a rather good description of the data 
(reduced x^ = 1.46) for both solutions reported in [29| . 
Including the new data in the database and slightly read- 
justing the parameters, we obtained an improved value 
of x^ = 1.26 without notable changes in the description 
of other photoproduction data. 

The highest spin state which was observed in photo- 
production reactions is the A^(2190)G'i7 resonance. This 
is a four-star resonance according to the PDG classifi- 
cation, and the behavior of the TriV elastic amplitude in 
this partial wave clearly shows the presence of a pole in 
the mass region 2100 - 2200 MeV/c^- The recent CLAS 
data on the reaction 7p — > Kh. [36| demand a contribu- 
tion from this state for a good description of the recoil 
asymmetry data. Thus, we expect to observe this state 
in the pion photoproduction reactions. Indeed, we found 
that the coupling of this state into the ttTV channel (which 
is fixed by the elastic data) is much stronger than that 
into the KK channel, and in the combined analysis of 
our database, we found a small but not negligible contri- 
bution from this state to the pion photoproduction cross 
section. The pole position of the G17 state was found to 
be 2170 — 1 155 MeV, and the mass scan is shown in Fig. [5] 
with a solid line. A minimum, fully compatible with that 
observed in the CLAS 7p — > ifA data is clearly seen. 

One of the new states which was observed in the earlier 
analysis of the pion and 77 photoproduction data [37[ was 
a D13 state with a mass of 2I66I5Q MeV/c^ and a width 
of 300 ± 65 MeV/c^. In the combined analysis of the 
present database, this state was optimized with a mass 
of 2130 ± 20 MeV/c^ and a width of 310 ± 30 MeV/c^. 
The mass scan of this state for the presented pion pho- 
toproduction data is shown in Fig. [9] as a dashed line. It 
is seen that the minimum corresponds well to the mass 
of the state found in the combined analysis, but the data 
are more tolerant to a lower resonance mass. 
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FIG. 9: Mass scans of the Gn state (solid line), the second 
_Di5 state (dotted line), and the fourth D13 state (dashed line), 
which show the deviations in fit quality from the optimum fit 
(Ax^) for a series of mass values. 



For the first observation of the N{2Q7Q)Di5 state in the 
analysis of the jp -^ p-q differential cross section [331 > i^ 
was shown that this state can also contribute to the pion 
photoproduction cross section. In the present analysis, 
the TriV coupling of this state was found to be similar to 
the rjN coupling. The mass scan, shown in Fig. [S] with 
a dotted line, produced the minimum slightly above the 
mass value defined by the combined fit. However, the 
data can accept well the mass of the state between 2060 
and 2120 MeV/c^. 

Higher-mass F-wave states can also contribute to pho- 
toproduction and were included as multichannel relativis- 
tic Breit-Wigner amplitudes in the PWA ^M,- We have 
performed mass scans of the 7V(2000)i^i5, iV(1990)Fi7, 
and Af(1890)S'ii states, which were observed in the anal- 
ysis of the 7p — > A'A CLAS data [30|. Candidates for 
these F-wave states are listed as two-star resonances by 
the PDG; a third ^n resonance has been assigned a one- 
star status [23]. Evidence for a higher- mass Fis-state 
around 2000 MeV/c^ was also recently reported by the 
CLAS collaboration in cj photoproduction off the pro- 
ton |38| . However, no minima were observed in our mass 
scans for the description of the present data. Further 
single- and double-polarization observables are needed to 
constrain the scattering amplitude and to extract weakly 
contributing resonances. 

In a final comment, we note that our description of the 
data and the current SAID solution (SP09) [30| underes- 
timate the most-forward data at all three photon energies 
presented in Fig.[Sl The overall systematic error can cer- 
tainly account for some of the observed discrepancy, but 
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some physics may be missing. This needs to be further 
addressed in future analyses. 



N{2080)Di3, and iV(2070)Di5 high-mass resonances in 
its description of the data. 



VIII. SUMMARY 
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In summary, we have presented data on the photo- 
produced TT*^ cross section in the reaction jp -^ pir^. 
The continuous beam from the ELSA accelerator and 
the fiber detector of the tagging system provided tagged 
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results are in good overall agreement with previous mea- 
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the models at very forward angles for large photon ener- 
gies. The Bonn-Gatchina PWA found the iV(2190)Gi7, 
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